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Abstract 

An iterative computer method is described for 
identifying boiler transfer functions using fre- 
quency response data. An objective penalized per- 
formance measure and a nonlinear minimization tech- 
nique are used to cause the locus of points genera- 
ted by a transfer function to resemble the locus of 
points obtained from frequency response measure- 
ments. Different transfer functions can be tried 
until a satisfactory empirical transfer function of 
the system is found. To ill ustrate the method, some 
examples and some results from a study of a set of 
data consisting of measurements of the inlet imped- 
ance of a single tube forced flow boiler with in- 
serts are given. 

Introduction 

Studies of the dynamic stability of single tube 
boilers have been conducted at the NASA Lewis Re- 
search Center. The aim of these studies was to es- 
tablish dynamic system stability criteria. The ex- 
perimental part of the program consisted of measur- 
ing the frequency response of the system. (1, 2, 3) 

As part of the analytical program a method was de- 
veloped to identify transfer functions from these 
measurements. This report presents the details of 
this method. Typical results obtained by applying 
the method to the boiler frequency response data 
are also presented. 

The frequency response data consist of the am- 
plitude ratio and the phase angle difference of two 
parameters at various frequencies. The frequency 
response is considered to be identified if a trans- 
fer function that fits the frequency response data 
is obtained. Various approaches can be used to 
analyze these data to obtain a satisfactory trans- 
fer function. 

If nothing were known about the dynamic system, 
a heavily empirical approach could be adopted. For 
example, it could be assumed that the system trans- 
fer function is composed of pure time delays e - ^, 
gains ki, and associated time constants Ti- It 
could further be assumed that the transfer function 
Z had the general form 



i=l 


where s = jo>. This curve fitting approach does 
not have generality. The constants are not related 
by theory to the steady-state condition and must 
be calculated for each test measurement. 

Another approach to the problem would be to 
attempt a rigorous solution of the partial differ- 
ential equations governing the boiler. This ap- 
proach assumes all important processes are known. 


The approach taken herein was to develop a 
method that could yield an analytical model. The 
method is used to obtain a transfer function with 
the following characteristics: The transfer func- 

tion should describe the measured data; the trans- 
fer function should preferably be derivable from a 
physical theory; end the dime nsionless gains, time 
delays, and time constants in the transfer function 
should be derivable from steady-state measurements 
or physical properties. This report is concerned 
with the development of the method and not with 
relating results to fundamental physical principles. 

In the next section the general method is ex- 
plained. In the then following section, the results 
of applying the procedure to boiler inlet impedance 
measurements are given. The last section contains 
a further discussion of this work. 

Method 

The success of the method is determined par- 
tially by the extent the locus of points generated 
by the transfer function (LFG) resembles the locus 
of points obtained from frequency response measure- 
ments (LTM) . The locus of points in either case 
is obtained by connecting the frequency response 
points in the order of increasing frequency. While 
it is necessary that this resemblance criterion be 
satisfied, resemblance is not a directly measurable 
quantity. A secondary criterion called an objective 
performance measure is used to provide a quantity 
that is directly measurable. 

A description of the manner in which this per- 
formance measure is used to investigate resemblance 
between LIM and LFG follows. Each physical descrip- 
tion (based on experimental results or fundamental 
physical principles) thought to describe the funda- 
mental dynamic phenomena is used to produce a system 
transfer function corresponding to the one measured. 

Each physical description yields estimates of 
the gains, time delays, or time constants in the 
transfer function. It also yields information on 
the possible ranges for these constants. For each 
transfer function the constants are selected to 
minimize the positive, quadratic performance meas- 
ure. The form of the performance measure used also 
ensures that all the constants will be within thier 
proper limits. In this manner each analytical 
model can be compared to measured frequency response 
data. 

The objective performance measure indicates 
which physical description is most applicable. The 
best LFG can be checked for resemblance by plotting 
and comparing with the LFM. 

In this manner a transfer function is con- 
structed by an iterative procedure. At each step 
the results can suggest modifications in the analy- 
sis that are necessary for improvement. This pro- 
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cedure can be used to pinpoint the inadequacies in 
an analysis. 

The basic performance measure S is 
J 

s = Y D z(si5 x) ■ ^^i ] 2 

iTi 

where s^ = jo> i and the (x) are the transfer 
function constants. The performance measure is the 
sum of the square of the magnitude of the separation 
of the measured impedance Z oks and the calculated 
impedance Z at «n frequencies where measurements 
are made. Previous reports (^, 5, 6) have suggested 
that a determination of the approximate system 
transfer functions and parameters can be made by 
finding the transfer function constants which mini- 
mize the performance measure. Standard computer 
techniques, to be discussed later, are used to find 
the constants (x) which minimize the performance 
measure. 


where 

0 for x(k,m)j < d(k,m)j 

1 for x(k,m)j > d(k,m)j 

This penalty function has the following character- 
istics, If the estimated parameter x(k,m)j for 
the k tjl theory and m*'* 1 steady-state condition is 
less than the estimated constraints, the penalty p 
is zero. If the estimated parameter x(k,m)j is 
greater than the estimated constraint, the penalty 
is equal to the square of the difference between the 
penalty and constraint. The penalty due to these 
individual penalty functions was defined as the 
weighted sum of these penalty function. 

With this penalty function, the objective per- 
formance measure S p was defined to be 

N 

S (k,m) = S w (k,m) + w iPi (x) 
i=l 


The actual performance measure used herein was 
a modified form of the previously given basic defi- 
nition. First, the basic performance measure S 
was changed to a weighted performance measure S w . 
This was done for various reasons. The weighted 
performance measure provides a simple means to ob- 
tain a performance measure that applies over a limit - 
ed frequency range. As an example, if measurements 
were made over a range from 0.001 to 4.0 hertz and 
if only the range from 0.001 to 1 hertz was to be 
covered by the transfer function, the weighting 
factors for frequencies larger than 1 hertz could be 
set to zero. The weighted performance measure 
furthermore provides a simple way to minimize the 
effects of noise or of errors on the identification 
procedure. It also provides proper scaling. The 
weighted performance measure has the form 

S w (k,m) = ^ a(k,m) i pZ k (s i ,{x(m)}) - Z m ( Si )Q 
i=l 

The k index represents the k^ theory. The m 
index represents the m th set of data. The weight- 
ing factor is designated as a(k,m)^. 

Since the selection of values of the transfer 
function constants is automated, one further modifi- 
cation was made to the performance measure. This 
modification was made to insure that the constants 
selected in minimizing the performance measure are 
in the proper range. The estimates of the proper 
range for each constant are used to provide con- 
straints on the selection of parameters, and an ac- 
ceptable variation of parameters is obtained. Hence, 
a transfer function need not be discarded if the 
initial selection of parameters was incorrect. An 
adequate understanding of the process or a lack of 
knowledge of important steady-state parameters need 
not prevent the recognition of a suitable transfer 
function if the range of the parameters can be 
estimated. 

The selection of parameters was constrained by 
the penalty function 

[p z (x) = d(k,m)j - x(k,m)j] 2 6j 


The weighting factor Wi is selected to scale the 
penalty such that each penalty will have an 

appropriate effect on Sp. 

The problem of minimizing a nonlinear function 
like the performance measure can be solved by a 
search technique. The basic methods available are 
discussed in chapter 6 of Eef. 7. To provide maxi- 
mum flexibility in the choice of Z, a search tech- 
nique which does not require evaluation of deriva- 
tives was chosen. The search technique used was 
that of Powell. ( Q ) The computer program used was 
adapted from Bef. 9. 

When the penalized performance measure and the 
computer methods mentioned previously are combined, 
the computer method shown schematically in Fig. 1(a) 
results. This method is used to find out if a par- 
ticular choice of transfer function fits the meas- 
ured frequency response data. 

The procedure used to obtain a satisfactory 
transfer function is shown in Fig. 1(b). Only the 
step in which Sp(k,m) is calculated is programmed 
on the computer. For the first transfer function it 
is best to select one based on a simple first-order 
theory that yields an approximation to the LCH4 and 
that uses the steady-state data available to predict 
the transfer function constants and constraints. At 
this point it is best to select the data that are 
easiest to work with and to limit the frequency 
range over which agreement is expected. This ini- 
tial transfer is denoted as Zp. The set of param- 
eters obtained is denoted as (x(l,m)) and the set 
of constants is (d(l,m)). The minimum value of 
Sp(l,m) can then be found. It should be noted that 
in this expression S p may be minimized to yield a 
transfer function Zq which is optimum over a 
limited frequency range. 

For the following iterative steps it is neces- 
sary to examine previous results to locate terms 
that would extend the frequency range, terms that 
would enable the theory to be extended to more 
steady-state operating conditions, and terms that 
were not used because of insufficient steady-state 
knowledge or physical knowledge. After the func- 
tional form of these neglected terms is found, the 
value of the parameters in these terms and the con- 
straints on these parameters are then estimated. 
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These terms when combined together produce various 
theoretical transfer functions. The k th theoret- 
ical description yields transfer function Z^. Also 
from the k™ theoretical description and for each 
steady-state condition, a set of constants (x(k,m)} 
and a set of constraints (d(k,m)} are available. 
The minimum value of Sp(k,m) is found. The re- 
sulting set of optimized transfer functions is 

compared. 

The results of each change in the theoretical 
transfer function, change in frequency range, or 
change in constraints on the previous theoretical 
transfer function are readily available using the 
conputer method. 

The procedure described previously was used to 
obtain a transfer function for the inlet impedance 
of a single tube, forced flow boiler with inserts 
from inlet impedance measurements. The following 
section discusses some of the results of this work. 
The aspects discussed in the next section are those 
not concerned with relating the transfer function 
obtained to a set of fundamental physical princi- 
ples. 


Application 


and (3) putting the equations in dime nsionless form. 
The di m ensionless form is obtained by dividing by 
scale factors P and W. The scale factors are 
not defined individually since in the final equation 
they appear as W/P. The quantity P/W is defined 
to he 3.28x10^ newtons -seconds per square meter - 
kilogram. The result is a set of equations connect- 
ing perturbations in parameters. The coefficients 
of the parameters are the subsystem transfer func- 
tions. Each subsystem transfer function represents 
a physical process. 


The perturbation parameters chosen to describe 
the inlet impedance are the following: the varia- 

tion of inlet pressure AP.^, the variation of inlet 
mass flow rate AW^ n , the variation of vapor pres- 
sure beyond the subcooled region APl SCj the vari- 
ation of subcooled length £L SC , and the variation 
in vapor mass flow rate beyond the subcooled region 
*WL S c- 


These five perturbations are connected by the 
following four equations: 



E(s)- 


ZM,- 




(1) 


f 



This section initi all y indicates how the contri- 
butions of various transfer functions (each attrib- 
uted to a different physical process) are separated 
in the formulation of the system transfer function. 
This is done so that each can be considered inde- 
pendently. The transfer function for one process 
can then be changed until the system transfer func- 
tion is improved while maintaining the transfer 
functions for other processes constant. Next, some 
typical results obtained by using the performance 
measure, computer techniques, and the iterative 
procedure described in the previous section are 
given. Last, the transfer functions which were 
found by this method are discussed. 

The analyzed data(2) were obtained experimen- 
tally by driving sinusoidally the open area of a 
valve in the feed system about a mean area and 
measuring the pressure and flow at the boiler inlet. 
The pressure and flow signals were analyzed by a 
frequency-response analyzer. This analyzer com- 
puted the magnitude and phase of the sinusoidal con- 
tent of the pressure and flow signals at the driving 
frequency. The magnitude of the boiler inlet im- 
pedance at each frequency is equal to the amplitude 
ratio of the inlet pressure to inlet flow signals. 
The phase of the boiler inlet impedance at each 
frequency is equal to the difference between the 
measured phase angle of inlet pressure and inlet 
flow. 

The theoretical model used to describe the 
measured boiler input impedance data is a function 
of a combination of physical processes. ' The dif- 
ferent physical processes are assumed to be tied 
together by a set of equations. While the detailed 
models of each physical process are not known ade- 
quately, the set of equations which forms the frame- 
work that ties the models together is assumed to 
be known. 

The equations and subsystem transfer functions 
are derived from the governing equations for the 
boiler by (l) simplifying and linearizing for small 
perturbations about a known steady-state operating 
condition, (2) Laplace transforming the equations. 


^sc . /»in , 

— — = G(s)-c^- + H(s)— =5 
L sc W P 


<2W t 


£W. 


AL„ 


AP T 


d 2 s 


- F(s)- 


AP, 


L sc 


(2) 

(3) 

(4) 


Equation (l) relates the perturbation in pressure 
drop across the subcooled region to the inlet mass 
flew rate through an inertial term E(s). Equa- 
tion (2) states the following: 

1. Changes in subcooled length (for constant 
vapor pressure) are due to perturbations in inlet 
mass flew rate effecting heat transfer (the trans- 
fer function for this physical process is repre- 
sented by the term G(s)). 

2. For constant inlet mass flow, perturbations 
in vapor pressure effect the subcooled legnth (the 
transfer function for this physical process is rep- 
resented by the term H(s)). 

3. These two effects can be added for small 
perturbations. Equation (3) states that changes in 
mass flow rate can be related to mass storage 
through changes in subcooled length (represented by 
the term 02 s ) and in pressure (represented by the 
term F(s)). Equation (4) relates variations in 
exit flow and exit pressure by the constant R. 

The inlet impedance derived using these equa- 
tions is 

AP in /P R(1 - e 2 sG(s)) 

Z(s) = = T= =r + E(s) 

£W in /W L 1 + R ( 0 2 sH ( s ) + F ( S ))J 

( 5 ) 

The basic block diagram for the system is shown in 
Fig. (2). 
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Each subsystem transfer function and thus each 
separate physical process can be investigated inde- 
pendently by proceeding in the previously stated 
manner. Thus, different theoretical models for 
each physical process can be checked. Also, differ- 
ent mathematical forms of the subsystem transfer 
functions can be tested. These forms can be those 
determined from physical theory or from transfer 
function plots. The method of this report was 
then used to determine the subsystem transfer func- 
tions G(s), E(s), H(s), and F(s). In the follow- 
ing section the results of a study of data taken at 
a condition of high vapor exit quality are dis- 
cussed. 


High Quality Vapor Exit Condition 


The time constant 6 g is designated as constant 
X}_. The time constant x 4 corresponds to the 
inertia of the subcooled fluid. It is given in 
the dimensionless form used here as 


x _ fisc W 
4 A P 


( 8 ) 


The gain Xg corresponds to the change in subcooled 
length with a change in vapor pressure. In nondi- 
men sional form it is 


Xg = 




SPl sc Lsc 


0) 


The case discussed is run 4 from Eef. 2. The 
data were obtained at a steady-state condition of 
high vapor exit quality (99 percent). The boiler 
inlet impedance data are plotted in Figs. 3(a) and 
(b). Figure 3(a) contains the data measured over 
the frequency range 0.04 and 1.0 hertz, and 
Fig. 3(b) contains plots of the data over the fre- 
quency range 0.25 to 4.0 hertz. The basic charac- 
teristic form is that of a spiral. Over the fre- 
quency range below 1 hertz (fig. 3(a)) the ampli- 
tude of the transfer function decreases as the 
phase angle rotates clockwise. Above 1 hertz 
(fig. 3(b)) three important changes in the form 
occur. Between 1.6 and 2.5 hertz a small loop 
occurs. Beyond 2.5 hertz after the small loop the 
following loop is larger. Also above 4 hertz the 
data are not symmetric about (0,0). Instead, the 
data are symmetric about a line through the real 
axis. Also, the spiral seems to be pulled up along 
a line parallel to the imaginary axis. Several 
theories are now compared with the data. 

Theory 1 . - The initial transfer function used 
to provide the starting point in the iterative 
model building process was the theoretical model 
described in Ref. 11. This model is designated as 
theory 1. This theory, which is based on many sim- 
plifying assumptions, provides an adequate estimate 
of the boiler impedance over a limited frequency 
range for those steady-state conditions where the 
theoretical simplifications agree with the actual 
steady-state conditions. The theory was not meant 
to be applicable to data taken at a condition of 
very low exit quality or very high exit quality. 


The model described in Ref. 11 used the follow- 
ing functional forms for the subsystem transfer 
functions G(s), E(s), H(s), and F(s): 

E(s) = x 4 s (6a) 


-0 2 s 


, 1 - e 

(6b) 

e 2 s 

H(s) = x 2 

(6c) 

F(s) = 0 

(6d) 


The time constant 9^ is the subcooled dead time 
and is given by 


L scP A 


The boiler resistance previously denoted R is 
henceforth defined as a constant: 

x 3 = R (10) 

To obtain a form of the inlet impedance equation 
that is more clearly related to more cctnmon forms, 
the following transformation of Eq. (5) was made. 
Using Eq. (10) and defining 

QUA = 1 - x 1 sG(s) (11) 


give 


X3QUA 

1 + x 3 x 1 sH(s) + x 3 F(s) 


E(s) 


( 12 ) 


The first transfer function used Ln the pro- 
cedure diagrammed in Fig. 1(b) was derived from 
Eqs. (11) and (12) and Eqs. (6a) to (d). The re- 
sult is 


Z = 


x 3 e 


-XqS 


1 + 


X 3 X 2 X 1 S 


Note that from Eq. (11) and (6d) 

-X-, s 
QUA-l = e 1 


(13) 


(14) 


Coefficients for this model were selected using 
the method described in the first section. The re- 
sulting function is plotted for frequencies less 
than 0.3 hertz in Fig. 3(c) and for frequencies 
greater than 0.3 hertz in Fig. 3(d). Examination 
of Fig. 3(c) shows that this model does describe 
the data in the range below 1 hertz. The data and 
the theory 1 locus both resemble a spiral. The 
amplitude of the transfer function decreases as the 
phase angle rotates clockwise. Examination of 
Fig. 3(d) shows that theory 1 fails in the frequency 
range above 1 hertz. Ths sizes of the loops calcu- 
lated decrease uniformly while the measured transfer 
function shows a smal 1 loop followed by a large one. 
Also in the same figure the calculated spiral is 
seen to be symmetric about the imaginary axis where- 
as the measured transfer function is seen to be sym- 
metric about a line parallel to the imaginary axis 
on the positive real side. An aspect of the 
theory 1 curve and the data that are in agreement 
is the appearance of loops. 


W 


( 7 ) 


The relation is now shown between Eq. (13) and 
the theory 1 curve in Figs. 3(c) and (d). The spi- 



t> 






f. 


ral shown in Fig. 3(c) is generated by the term 


1 + XgXgXjS 


(15) 


The decrease in amplitude occurs becuase the magni- 
tude of 


glance at Fig. 3(f) indicates that the theory is 
not successful since the resulting curve does not 
resemble the data at high frequencies. The theory 2 
loops decrease in size rapidly as the frequency in- 
creases. Also, no small loop followed by a large 
one occurs. This modification is an example of how 
the extension along well known lines of an idea may 
result in an unsuccessful transfer function. The 
final version is discussed next. 


x 3 

1 + x 3 X2X lS 

decreases with frequency. 


(16) 


The change in phase angle with frequency is 
seen as a rotation in Figs. 3(c) and (d). The 
rotation is due to a gentle change in phase angle 
from Eq. (16) combined with a rapid change in phase 
angle due to the term 


e‘ X l S (17) 

The addition of the imaginary term 

x 4 s (18) 

to term (15) causes the spiral to be pulled up along 
the imaginary axis. This effect produces the loops 
seen in the theory 1 curve of Fig. 3(d). Theory 1 
can be modified. The next sections discuss one 
type of modification. 

Theory 2 . - As an example of one step in the 
iteration procedure, the results obtained by modi- 
fying the previous model to include the effect of 
wall storage on the heat transfer process are dis- 
cussed. The wall storage modification is typical of 
a modification that can usually be made frcm a more 
detailed analysis of the physical processes in a 
system. 

The modification results in a change in the 
definition of QUA. This modification will be called 
theory 2. The modified term was 

QUAg = pe -?v + 1 - M. (19) 


where 


Theory 5 . - The model found to represent the 
high vapor exit quality data best consisted of the 
following set of subsystem transfer functions: 


-x n s 

QUA = e ~ 

(23) 

s 

V S + — ; 

J- 1/X 1Q + S/Xg 

(24) 

“XcS 

H = x^se 0 

(25) 

F _ V 
2 2 

(26) 

s +■ XgS x 7 



The resulting transfer function, which is 
plotted in Figs. 3(g) and (h), is designated as 
theory 3. This model has all the important charac- 
teristics of the data at all frequencies. The 
transfer function has the basic spiral shape below 
1 hertz (fig. 3(g)). It also has the small loop 
between 1.6 and 2.5 hertz. And the larger loop 
between 2.5 and 4 hertz also exists (fig. 3(h)). 
Furthermore, the high frequency data are symmetric 
about a line through the real axis. The inertia 
effect that pulls the data up parallel to the imag- 
nary axis is also apparent. . 


The transfer function is 


Z 1B 


-XcS 

1 + x-^XgXgse ° + 


X 3 X 8 S 


+ x 4 s + 


s*- + 


XgS + x^ 


x 1 s (l ' 

(20) 

1 + x 12 / x ll + X 1 S / X H 

, x l s + x 12 * x 12 

(21) 


1 + X 1 s/x^ J 


Thus, using Eqs. (12), (6a), (6d), and (19) yields 


(pe ^ + 1 - p)x 3 


^lA 


1 + 


X 1 X 2 X 3 S 


( 22 ) 


s 

l/ x 10 + s / x 9 


(27) 


The appearance of the s m all loop before the 
large one (fig. 3(d)) is due to the additional term 
F(s). The magnitude of the term 


1 + xj_xgx 3 se 


-V + 


x 3 x 8 e 


+ XgS + 


becomes very small at the frequency vhe 


(28) 


A range of parameter values for x^, x, g, and 
x,, was selected. Also, values of the coefficients 
for the specific operating condition of the example 
was chosen. The resulting transfer function is 
shown in Figs. 3(e) and (f). 

Examimation of Fig. 3(e) shows that the 
theory 2 curve describes the data at low frequencies 
as well as the theory 1 curve. The procedure used 
has minimized the performance index. Still, a 


This decrease in magnitude causes the loop that 
occurs between 1.6 and 2.5 hertz to be small. The 
high frequency data are symmetric about a line 
through the real axis parallel to the imaginary 
axis due to the term 
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Appendix - Symbols 


1 / x 10 + s / x 9 


(29) 

A 


cross-sectional area of toiler tube, 


in E(s). At high frequencies including the term d 

causes E(s) to become 


constraint on value of transfer function 
constant 


E(S) = X 4 S + Xg 


(30) E(s) ■ 


subsystem transfer function 


In conclusion, it should be noted that before F(s) 
Eq. (27) was determined no equation existed which 
gave good agreement with the form of the high vapor f 
exit quality data for the frequency range above 
1 hertz. G(s) 


subsystem transfer function 
frequency, Hz 

subsystem transfer function 


The three theoretical models used in the study 
of the high vapor exit quality example are summar- 
ized in table 1. The set of coefficients that were 
obtained using the method with each model and the 
data are given in table 2. Some additional results 
of the study are presented next. 

Low Quality Vapor Exit Condition 

•Another result of this study was the discovery 
that a modification of one term in Eq. (27) would 
enable that equation to be used to interpret the 
inlet impedance data taken at a condition of low 
vapor exit quality. The modification was to use the 
following form of F(s) instead of Eq. (26). 


H(s) 

I 

j 

k 

Lgc 

^scAsc 

p 

p 


subsystem transfer function 
inertia, sec 

gain 

subcooled length, m 
perturbation in subcooled length 
pressure, N/m 2 
mean pressure, N/m 2 


F(s) 


s x 3 

1 + s x 1 1 


AP/P perturbation in pressure 

p penalty function 


The excellent agreement possible with this B 

transfer function is demonstrated for a typical 
low vapor exit quality (20 percent) case( 2 ' (case 1) S 
in Fig. 4 where a comparison of the measured and 
calculated boiler inlet impedance is made. S p 


subsystem transfer function 
performance measure 
penalized performance measure 


The procedure described in the report enabled 
a determination of the transfer function form to be 
made for low exit vapor quality and high vapor exit 
quality. Thus it has proved to be valuable in ex- 
tending and increasing the understanding of the dy- 
namics of boilers. 

Concluding Remarks 

The search technique used permitted the param- 
eters which were derived from the most certain 
theoretical aspects of the problem to stay relative- 
ly unchanged during the search for a minimum of S p . 
A good example was the term xg which is calculated 
from the slope of the pressure drop against flow 
curve. (H) This term determined the magnitude of 
the transfer function at low frequencies, and it 
varies by less than 2 percent in a typical computer 
calculation. 

The parameters which were less certain took on 
a larger range of values during a search. The dif- 
ficulty is that the Powell method generally provides 
only a local minimum and this means that the solu- 
tion is not inherently unique. The theoretical re- 
sults are still considered valid since they apply 
to a range of experiments with different steady- 
state conditions and boiler geometries.' The value 
of the constants obtained are also physically 
reasonable since the use of penalty functions pro- 
vides estimates which are possible. 


s 

T 


W 

w 

AW/W 

w 

(x) 

Z 

Z = (P/w) 



weighted performance measure 
Laplace operator, sec”'*’ 
time delay, sec 
mass flow, kg/sec 
mean mass flow, kg/sec 
perturbation in mass flow 
weighting factor in penalty function 
set of transfer function constants 
transfer function 

impedance value relative to which all 
transfer functions are compared, 
scale factor 3.28x10 s (n-sec/m 2 kg) 

weighting factor in performance measure 

weighting factor for penalty 


P 


J 


% 
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TABLE 1. - COMPARISON OF FUNCTIONAL FORMS 
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TABLE 2. . DIMENSIONLESS TRANSFER FUNCTION 
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Figure 1. - Flowcharts illustrating identification procedures. 
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(b) Boiler data for frequencies above 0. 25 hertz. 












(d) Comparison of data and theory 1 for frequencies above 3 hertz. 
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(f) Comparison of data and theory 2 for frequencies above 0.3 hertz. 
Figure 3. - Continued. 









